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ABSTRACT 

We use two-dimensional axisymmetric magnetohydrodynamic simulations to compute steady-state 
solutions for solar-like stellar winds from rotating stars with dipolar magnetic fields. Our parameter 
study includes 50 simulations covering a wide range of relative magnetic field strengths and rotation 
rates, extending from the slow- and approaching the fast-magnetic-rotator regimes. Using the simula- 
tions to compute the angular momentum loss, we derive a semi-analytic formulation for the external 
torque on the star that fits all of the simulations to a precision of a few percents. This formula pro- 
vides a simple method for computing the magnetic braking of sun-like stars due to magnetized stellar 
winds, which properly includes the dependence on the strength of the magnetic field, mass loss rate, 
stellar radius, suface gravity, and spin rate and which is valid for both slow and fast rotators. 
Subject headings: magnetohydrodynamics — stars: magnetic field — stars: rotation — stars: solar- 
type — stars: winds, outflows 



1. INTRODUCTION 

The evolution of the rotational properties of sun-like 
stars exhibit numerous trends that still lack accepted, 
quantitative explanation. For example, during the pre- 
main-sequence phase, there is not yet a comprehen- 
sive picture that explains the observed wide range of 
rotation rates, nor the fact that a large fract i on of 
stars rotate relatively slowly (e.g.. iHerbst et al.l 120071: 
llrwin fc Borrvien 1200a iMatt fc Pudritzl l2005bllaf r Dur- 
ing the main sequence phase, we are still trying to un- 
derstand, for example: the structure and existence of 
multiple "sequences" apparent in the rotation- period- 
versus - mass diagrams of young star clusters (e.g..lBarnesl 
20035 Ihwm fc Bouvierl f2009; Mcibom et al. ||201J; the 
apparent "sat uration" of angular momen tum loss in 
fast rotators (jMacGregor fc Brennerl *~ 
all, se cular spin- down of stars (jKraftl 



1991); the over 



1967t iSkumanichl 



[19721: ISoderbloml Il983h : the possibility of using rota- 
tional properties to measure stellar ages ( "gyrochronol- 
ogy" :lBarnesll2003al l2010h iMamaiek fe Hiilenbrandl l2008l: 
iMeibom et al.ll2011l;lEpstein fe Pinsonneaultll2012f ): and 
the correlation of stellar activity with rotation, as well 
as th e "saturation" of this act i vity in fast rotators 
(e.g.. iSaar fc Brandenburg! 11999 1 : iPizzolato et all 120031 : 
iReiners et al.ll2009t iWright et al.ll201lD . 

We know that magnetized stellar winds are impor- 
tant for extracting ang ular momentum from stars dur- 
ing the main sequenc e (iParkerl 119581 : iSchatzmanl Il962t 
I Weber fc Davisl 119671 : iMestell 119681) and likely dur- 
ing pre-main-sequenc e fHartmann fc MacGregor! 11982b 
iHartmann fc Staufferl 119891 : IMatt fc Pudritzl l2005al) . 
Thus, a prescription for calculating the stellar wind 
torque as a function of stellar parameters is a cru- 
cial ingredient in models for the ro t ational evolu- 
tion of stars (e.g . . iBouvier et al.l 119971: [Bouvierl 120081 : 
iDenissenkov et al.ll2oTob IMatt et al.ll2012t ). 

Reliably computing the stellar wind torque requires 



knowledge of the wind acceleration profile and the mag- 
netic field ge ometry above the surface of the star (e.g., 
IMestell 11984). Until a few years ago, the only formula- 
tions available for computing stellar wind torques (e.g., 
Kawalcr 1988) were based upon analytic or semi-analytic 
calculations that necessarily relied upon several simplify- 
ing assumptions, such as that of spherical symmetry and 
the a priori specification of the magnetic geometry, flow 
acceleration profile, or both. However, in real winds, all 
of the flow properties are determined by a complex inter- 
action between the wind driving physics, stellar rotation, 
and dynamical interaction between the wind and mag- 
netic field, and these significantly deviate from spherical 
symmetry. Multi-dimensional, magnetohydrodynamical 
(MHD) simulations provide a reliable method for com- 
puting dynamically self-consistent wind solutions, but 
large parameter studies are needed, in order to quantify 
how the physics scale with various parameters. Thus, al- 
though there exists a large body of literature on the sub- 
ject, reliably calculating stellar wind torques for a range 
of stellar and wind properties, and in a way that is useful 
for stellar evolution calculations, remains a challenging 
problem. 

Using simula ti ons of simplified, solar-like winds, 
IMatt fc Pudritzl (|2008l hereafter MP08) carried out a 
small parameter study to determine the dependence of 
the wind torque on the strength of the magnetic field and 
mass outflow rate in the winqj. Their formulation for the 
torque (equation 15 of MP08) differed significantly from 
previous formulations, notably in the value of the expo- 
nent in the power-law dependence of the torque on stellar 
properties, such as magnetic field strengh, mass loss rate, 
and stellar radius. This study adopted many of the same 
assumptions as previous analytic work — ideal MHD, a 

1 Although the results were presented in the context of pre-main- 
sequence angular momentum loss, the simulations of MP08 are also 
appropriate for main sequence stars. 
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rotation-axis-aligned dipolar magnetic field, solid body 
stellar rotation, and spherically symmetric thermody- 
namic properties at the stellar surface — but the simula- 
tions did not require any assumptions about the kinemat- 
ics of the flow nor how the magnetic geometry was mod- 
ified by the flow. Thus, the MP08 torque is the most dy- 
namically self-consistent formulation for the torque from 
sun-like stars to date, and the implications for stellar 
evolution are still being explored. 

At the same time, this formulation is derived from 
simulations with variations only in the magnetic field 
strength (relative to the mass loss rate and surface grav- 
ity). It does not fully capture the effects of different rota- 
tion rates, different thermodynamic (or energetic) prop- 
erties of the wind, nor more complex magnetic geome- 
tries. Given the importance of computing stellar wind 
torques for a range of stellar ages (and thus a range of 
rotation rates), the natural next step is to extend the 
parameter study of MP08 to include variations in both 
magnetic field strength and stellar rotation rate. In this 
letter, we present such a parameter study and, from these 
results, derive the most generally applicable stellar wind 
torque formula to date. 

2. SIMULATION METHOD AND PARAMETER STUDY 

We use the simulation code and method described in 
MP08 to compute steady-state wind solutions for sun- 
like stars with a dipole magnetic field. We briefly de- 
scribe the method here, and the reader will find further 
details in MP08 and references therein. The code solves 
the equations of ideal MHD under the assumption of ax- 
isymmetry and a polytropic equation of state (P oc p 7 ). 
In each simulation, the numerical grid is initialized with a 
spheric ally symmetr ic, thermally-driven Parker wind so- 
lution ()Parkerlll958l ). plus analytic dipole magnetic field. 
Once the simulations begin, the wind solution relaxes to 
a steady state resulting from a dynamical balance be- 
tween the accelerating wind and rotating magnetic field. 
The steady-state solution is entirely determined by the 
conditions that are present at the base of the wind (the 
"surface" of the star). 

For a given initial magnetic geometry, unique wind 
solutions are determined by dimensionless parameters, 
which can be given as three velocity ratios — va/v csc , 
Cs/^csc, and / specified at the surface and equator of the 
star — plus the adiabatic index 7. Here, va is the Alfven 
speed, Ucsc the gravitational escape speed, c s the ther- 
mal sound speed, and / is the equatorial rotation speed 
divided by the breakup speed, 

v A /v csc = B r (47rp r )- 1 / 2 (2GM r /R r )- 1 / 2 , (1) 

/ = nM /2 (GM^ 2 , (2) 

where £?» is the magnetic field strength at the stellar 
equator, p» the mass density at the base of the wind, 
G Newton's gravitational constant, M* the stellar mass, 
.R* the stellar radius, and the (solid body) angular 
rotation rate of the stellar surface (— 2n/P* : where P* is 
the rotation period). 

The torque formulation of MP08 is based upon 9 sim- 
ulations with variations in the parameter va/v csc . They 
also presented 5 simulations with variations in the other 
parameters, /, c s /v csc , or 7, which demonstrated that 
these parameters affect the torque in a way that is not 



Table 1 

Simulation Parameters and Results 



Case 


/ 


VA/Vesc 


T 


r A /R* 


1 


0.0000995 


0.0753 


27.5 


5.15 


2 


n rinnnrinc 

0.0000995 


0.301 


446 


9.42 


3 


0.0000995 


1.51 


13500 


19.5 


1 


U.UUU997 


0.0301 


4.11 


3.31 


5 


r\ nnnoov 

u.uuuyy i 


0.0753 


25.8 


5.04 


6 


0.000997 


0.301 


437 


9.58 


7 


r\ nt~\r\r\r\ r 7 

U.UUU997 


1.51 


13000 


19.4 


8 


U.UUo9o 


0.0753 


25.8 


5.03 


9° 


U.UU00O 


0.209 


215 


8.36 


10 


0.00393 


0.301 


436 


9.57 


11 


n nrno£' 

U.UUood 


0.418 


837 


10.9 


12 


n nnQOQ 
U.UUo9o 


0.953 


4900 


15.8 


13 


U.UUo9o 


1.51 


13000 


19.4 


14 


0.0101 


0.0753 


25.6 


5.01 


15 


0.0101 


0.301 


432 


9.50 


16 


0.0101 


1.51 


12700 


19.4 


17 


n nono 

U.UzUz 


0.0753 


24.9 


4.95 


18 


n nono 
U.UzUz 


0.301 


417 


9.31 


19 


n nono 
U.UzUz 


1.51 


11900 


19.2 


20 


n nono 

0.0299 


0.0753 


24.0 


4.86 


21 


nono 

U.Uz99 


0.301 


395 


9.04 


2J 


n nonn 
U.Uzy9 


1.51 


10900 


18.7 


23 


n a no 

0.0403 


0.0753 


22.8 


4.74 


24 


n a no 


0.301 


367 


8.73 


25 


n n a nQ 
U.U4Uo 


1.51 


9840 


18.1 


26 


n n a no 
U.U49o 


0.0753 


21.6 


4.62 


27 a 


n a no 

0.0493 


0.235 


212 


7.66 


28 


n a no 


0.301 


341 


8.44 


29 


n n a no 
U.U49o 


1.51 


8960 


17.5 


30 


0.0594 


0.0753 


20.1 


4.49 


31 


0.0594 


0.301 


312 


8.12 


32 


n n crn a 


1.51 


8020 


16.8 


33 


0.0986 


0.0301 


2.70 


2.82 


34 


nnor 1 

U.U9od 


0.0753 


14.6 


3.96 


35 a 


n nnoc 
U.U9oD 


0.209 


104 


5.98 


36 a 


n nno7 


0.254 


150 


6.44 


37 a 


nnor 1 

0.0986 


0.301 


211 


6.97 


38 a 


nnon 

0.0986 


0.363 


295 


7.53 


39° 


0.0990 


0.502 


580 


8.77 


40 a 


0.0986 


0.602 


857 


9.58 


41° 


0.0986 


0.940 


2130 


11.8 


42 a 


0.0986 


2.10 


10400 


16.8 


43 a 


0.0986 


3.01 


20500 


19.3 


44 


0.197 


0.0753 


5.39 


2.83 


45 


0.197 


0.301 


75.9 


4.93 


46 a 


0.197 


0.495 


212 


6.26 


17 


0.197 


1.51 


2120 


10.5 


48 


0.403 


0.753 


88.0 


4.61 


49 


0.403 


1.51 


370 


6.30 


50 


0.403 


3.01 


1470 


8.33 



a Same Case used in the study of MP08. 



captured by their fit formulation. Motivated by the fact 
that the study of MP08 is based on a relatively small 
number of simulations and by the importance of pre- 
cisely determining the torque for a range of stellar rota- 
tion rates, the present work extends the parameter study 
of MP08 to include 50 simulations sampling a large range 
in both UA/i>esc and /. For all simulations presented here, 
we adopt a dipole magnetic geometry and fix 7 = 1.05 
and c s /v csc = 0.222, as in the "fiducial" case of MP08. 

Table [1] lists the input parameters / and va/v csc for all 
50 cases of the present study. Twelve of the simulations, 
as indicated in the table, are identical to those used in 
MP08. 

3. SIMULATION RESULTS 
From the steady-state wind simulations, we compute 
the total mass loss rate (M w ) and the torque (r w ) di- 
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Figure 1. Parameter space explored. Each point represents a 
single stellar wind simulation. Symbols are as follows: asterisks 
represent cases with / < 0.02, circles are for / fa 0.03 and 0.4, 
triangles are for / fa 0.04, plusses are for / fa 0.05, squares are for 
/ fa 0.06, diamonds are for / fa 0.1, and X's are for / fa 0.2. 

rectly, by integrating the mass and angular momentum 
flux over a closed surface containing the star. Table [1] 
lists the value of 



2 n2 



T EE 



BiR 



(3) 



for each case. The quantity T is a dimensionless way 
of expressing the (inverse) mass loss rate resulting from 
the simulations, it is physically related to the simulation 
parameter (uaAw) > an d i t is similar to the "m a gnetic 
confinement paramet e r" of lUd-Doula fc Owockil (2002, 
lUd-Doula eFaLl 120081 . I2009D . In the following analysis, 
we treat T as an independent variable of the simula- 
tions, even though it is a result of the simulations. We 
prefer to work with T, as it contains the more observ- 
able/predictable quantity M w , as opposed to va/v csc , for 
which one must specify the density at the base of the 
corona p* (see eq. PQ ) . Furthermore, for a given value of 
psr, the mass loss rate depends sensitively upon the val- 
ues of c s /v csc and 7 (see MP08), as well as the particular 
choice of boundary conditions, which we hold fixed in all 
simulations. Thus, the consideration of T as a controlling 
parameter (instead of VA/v esc ) significantly mitigates the 
influence of c s /v esc , 7, and boundary conditions choices 
in the resulting torque formulation. 

Figure [T] illustrates the range of parameter space, in 
terms of / and T, spanned by all 50 simulations. Each 
symbol in the plot corresponds to a simulation listed in 
Table [TJ and the different symbol styles highlight differ- 
ent ranges of spin rates. The asterisks are for cases in 
the slow rotator regime, where the rotation has a neg- 
ligible influence on the wind dynamics. The remain- 
ing symbols represent cases in which rotation affects the 
speed and collimation, and consequently the efficiency 
of angular momentum loss, in the flow. For reference, 
the appropriate values for the solar wind are / « 0.004 
and T within a possible range of ~ I0 2-3 (e.g., for 
M w w 2 x I0 -14 M© yr _1 and an equatorial dipolc field 
strength of 1-5 Gauss). 

In order to express the resulting torques in a useful 
and general way, we consider the following. In a steady- 
state wind, under the assumptions of ideal MHD, the 
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Figure 2. Alfven radius (eq. [5]) versus T. The symbols have the 
same meaning as in Figure [l] Dotted lines show the fit function 
(eq. [6] and □) for / = 0.001, 0.1, 0.2, and 0.4, from top to bottom, 
which approximately correspond to a few of the rotation rates used 
in the parameter study. 

specific angular mo mentum extracted fro m the star is 
equal to fi*r A (e.g.. iWeber fc Davis! 119671 ). where ta is 
the "Alfven radius," the cylindrical radial location where 
the wind velocity equals the local Alvfen speed. In a 
three-dimensional flow, the net torque on the star can be 
written 

r w = M w ny A , (4) 

where r\ is an average value (|Washimi fc Shibatal[l993] ) . 
From the simulations, we compute r w and M w directly, 
so we use equation @ to define the dimensionless Alfven 
radius as 



R* 



1/2 



(5) 



which thus represents the square root of the dimension- 
less torque and is unique for a given set of dimensionless 
simulation parameters. Table [JJ lists the resulting values 
of ta/R* for each case. 

4. NEW TORQUE FORMULATION 

Figure [2] shows the value of the dimensionless Alfven 
radius versus T for the entire parameter study. The 
study of MP08 found the dimensionless Alfven radius 
was well- represented by a single power-law in T, for all 
cases with the same spin rate. In other words, 9 sim- 
ulations in MP08 were fit by ta/R* = KT m , where K 
and m were dimensionless fit parameters. In Figure [5J 
it is evident that the effect of rotation is to modify the 
value of K , but not to significantly affect the exponent m 
(represented by the slope in the log- log plot). Thus, we 
are able to fit the 50 simulations here, with variations in 
both T and /, with only one additional free parameter. 
Specifically, we find that all of the data are well fit by 
the function 



T A 



T 



(K$ + 0.5/ 2 )!/2 



(6) 



where Ki, K2, and m are dimensionless fit constants. 
The best fit values give: 



Kx = 1.30, 



K 2 = 0.0506, 



m 



0.2177 (7) 
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Figure 3. Dimensionless Alfven radius times T _m versus /. This 
Figure demonstrates the dependence of the Alfven radius on the 
stellar spin rate, for a fixed value of T. The symbols have the same 
meaningas Figure[l] and the dotted line shows the fit function (eq. 
[6] and 0). 

These values differ only by a few percent from the equiv- 
alent fit parameters of MP08 (m and K) but should be 
considered to be more precise here, due to the larger 
number of simulations fit. Also, the functional form of 
equation ([6]) simultaneously quantifies the effect of the 
wind magnetization (T) and the spin rate (/) on the ef- 
fective Alfven radius. The dotted lines in Figure [5] show 
equation ([6]) for four different values of / = 0.001, 0.1, 
0.2, and 0.4 and using the best fit values Each dot- 
ted line goes through the symbols that represent simula- 
tions with corresponding spin rates, illustrating how well 
equation ([6|) fits the simulation results. 

One can understand the functional form of equation 
([6]) as follows. The most important factor for determin- 
ing the Alfven radius is the strength of the magnetic field 
compared to the inertia in the flow, T. The MHD simula- 
tions self-consistently capture how the wind accelerates, 
how the magnetic field strength varies with distance from 
the star, and how much of the total surface magnetic flux 
will participate in the wind (the remaining flux exists as 
closed magnetic loops) . The fit values of K\ and m quan- 
tify how these processes depend upon the value of T, for 
a fixed rotation rate. 

For different rotation rates, the Alfven radius may be 
modified further because rotation can act as an addi- 
tional wind driving component. When the stellar rota- 
tion is very slow, the rotation has a negligible affect on 
the wind driving. However, for fast rotation, magneto- 
centrifugal effects increase the wind velocity. In order to 
quantify this effect, one can think of the wind speed as 
being proportional to a rotation-modified speed, 



K. 



2„,2 



! sc (^ 2 2 + 0.5/ 2 ). (8) 

In this sense, equation © is equivalent to the Alfven 
radius being a simple power law in equation ([3]), but 
with the factor of w csc being replaced by w mo d- The di- 
mensionless factor Ki determines at what spin rate the 
stellar rotation becomes dynamically important for the 
wind. 

The effect of the stellar spin rate on the Alfven radius 
(for a fixed value of T) is best illustrated by Figure |3l It 
is clear that when the spin rate is below a few percent of 



Figure 4. Dimensionless Alfven radius versus the quantities in 
square brackets in equation Jfj). This Figure demonstrates the 
dependence of the Alfven radius on T, for a fixed stellar spin rate. 
The symbols have the same meaning as Figure [l] and the dotted 
line shows the fit function (eq. [6] and 0). 

the breakup rate, the Alfven radius is independent of spin 
rate. Faster rotation decreases the Alfven radius. For the 
fastest spin rates in the parameter study (/ » 0.4), the 
Alfven radius is decreased by a factor of approximately 
30%, relative to the slowly rotating regime. This cor- 
responds to a factor of approximately 2 in the torque 
(r w cx r|). 

Figure, |4] shows the power-law dependence of the 
Alfven radius on T, for a fixed value of /. It is remark- 
able how well equation ([6]) fits the dimensionless Alfven 
radii determined from all 50 simulations. As can be seen 
in Figures HHH all data points lie within a few percent 
(the biggest outlier is off by 4%) of the fit function. This 
reflects the precision of the simulation method in deter- 
mining the values of r^l R* and T. 

By combining equations ((2]), ((3]), ©, and ((U, the 
torque on the star, due to the stellar wind, can be writ- 
ten 



r w = K(Bt m M 



2 £>4m ■\' / fl — 2mj^Am-\-2 



K 2 ■ f 

■"■1 g4mi|^l-2mp4m+l„,l-2m J_ 

V2 * w 



(2G) ? 



-Bz m M„ 



osc (K\ + 0.5/ 2 ) m 
Rl m+2 



2\m ' 



(9) 



where we have listed 3 different equivalent forms for con- 
venience. This formula is suitable for studies of the evo- 
lution of stellar angular momentum. 

5. DISCUSSION 

The torque formulation presented here differs signifi- 
cantly from analytic prescriptions that have be en widely 
used. Of particular note, the preferred model of lKawaien 
( 1988) results in a power-law dependence of the Alfven 
radius on T, with a power-law exponent of 0.5. In the 
present work, we find this exponent is approximately 
0.22, which results in a significantly different dependence 
of the torque on all stellar parameters. A comparable 
exponent of between 0.2 and 0.25 was found in s olar- 
like winds s i mulati ons of Washimi & Shibata ( 1993J and 
iPinto et all (|20ill) . and lUd-Doula et al.l (|2009i rfound 
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a similar exponent of 0.25 in simulations of radiation- 
driven winds from massive stars. The reasons and im- 
plications for m « 0.2 — 0.25, compared to the analytic 
work, was discussed in MP08. 

We have determined the effect of the stellar rota- 
tion rate on the Alfven radius. As evident in Fig- 
ure [21 there is a slow magnetic rotator regime, where 
the rotation rate does not i nfluence the Alfven radius 
(|Belcher fc MacGregorlll976l) . However, for faster rota- 
tion the increased wind acceleration by magnetocentrifu- 
gal effects acts to decrease the Alfven radius. This effect 
and transition is not captured in previous torque formu- 
lations. Note that our simulation parameter study does 
not extend fully into the fast magnetic rotator regime, 
where the magnetocentrifugal a cceleration completely 
domi nates the thermal driving ([Belcher fc MacGregorl 
Gill. Rather, even for the fastest rotation rates consid- 
ered here, the thermal wind driving was not negligible, 
as appropriate for solar-like winds. 

The torque formulation presented here provides the 
most physically realistic and precise calculation to date 
of solar-like angular momentum loss, and it is suitable 
for any studies of the evolution of angular momentum of 
sun-like stars. However, it has a number of limitations, 
which point the way for future work. First of all, the 
formulation is not valid in the limit of very weak mag- 
netic fields. In particular, when the size of the Alfven 
radius approaches the stellar radius, the angular momen- 
tum transport will begin to be dominated by other (e.g., 
viscous) effects, which are not properly included in our 
simulations. To take such effects approximately into ac- 
count, we suggest imposing a minimum value of r\> i?*, 
in situations where weak magnetic fields may be consid- 
ered. 

Secondly, the present study fixed the physical parame- 
ters that control the thermal driving of the wind, 7 and 
Cs/^osc (which are related to the coronal temperature and 
heating/cooling physics). MP08 included a few simula- 
tions that had variations in these wind driving param- 
eters, and we can also com pare the torque f ormu lation 
presented here with that of lUd-Doula et all (|2009f ) . de- 
rived for radiatively driven winds. Taken together, it is 
clear that a reasonable range in possible wind driving 
physics can affect the torque on the star by a factor of 
~ 2. It also stands to reason that the transition from slow 
magnetic rotator to the regime where stellar rotation is 
dynamically important in the wind (i.e., the value of Ki 
in eq. [6]) will depend upon the wind driving physics. 
It will be important in future work to be able to reli- 
ably predict how the wind driving physics systematically 
affect the torque. 

Finally, the present study assumed a rotation-axis- 
aligned, dipolar magnetic geometry at the stellar surface. 
This is justified by the fact that the largest-scale, global 
field has the most influence on the torque, but it is not 
yet clear how more complex (realistic) magnetic config- 
urations will change the scalings. MP08 presented two 
simulations with pure quadrupolar magnetic fields. For 
the same surface magnetic field strength B», the torque 
for the quadrupole case was reduced by an order of mag- 
nitude compared to the dipole case. They also noted 
that the power-law dependence of the Alfven radius on 
the strength of the magnetic field may be altered (their 2 
simulations suggested an exponent of 0.15, instead of 0.22 



for the dipole). Similarly, iPinto et all (|201lD simulated 
solar winds with a range of magnetic field complexity, 
noting a decreased Alfven radius with increased mag- 
netic complexity and suggesting a different power-law 
exponent when the field is significantly non-dipolar. It 
will be important in future work to systematically study 
how the torque is affected by complex, realistic magnetic 
geometries. 
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